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TIm coaputatlonal Mthods used to predict and optlain the therael*- 
structural behavior of aerospace vehicle structures are reviewed. In general* 
two classes of algorlthns* Inpllclt and explicit* are used In transient ther* 
■al analysis of structures. Each of these two anthods has Its own werlta. 

Eue to the different tlwe scales of the eechanlcal and thenal responses* the 
selection of a tine Integration aethod can be a difficult yet critical factor 
In the efficient solution of such problras. 

Therefore sized tine Integration eethods for transient thenal analysis 
of structures are being developed. The coeputer laplenentatlon aspects and 
minerlcal evaluation of these sized tlse lapllclt^ezpliclt algorlthas in 
ttersal analysis of structures are presented. A cosputatlonally-useful aethod 
of estlaatlng the critical tlse step for linear quadrilateral eleaent Is also 
given. Nuserlcal tests confln the stability criterion and accuracy 
characteristics of the sethods. The superiority of these sized tlse aethods 
to the fully Ispllclt aethod or the fully explicit aethod Is also 


deaonstrated 


Ovttr ctM l«at two d«cad«a» tlgnlf leant attantlon has baan davotad to tha 
davalopaant of lightaalght» durabla thanal protactlon sjataM (TPS) for 
futura apaca transportation syataas. fiasaarch prograaa ara eurrantly undar- 
way at tha Langlay Sasaarch Cantar to Invastigata various aatalllc TPS eoncapt 
[1]. Ona of tha proposad candidatas Is tha tltaniua ailtlwall tlla (saa [2] 
and rafarancas tharaln for a discussion). Early daslgn procaduras of tha TPS 
eoncapt Involvad both analytical and azparlwantal studies. In particular, a 
dag*ae of confidence has baan established In tha TPS concept due to tha daslgn 
studies by Jackson and Dixon [3] and Blair at al. [4]. 

A tltanltai aultlwall tile consists of alternating layers of suparplas** 
tlcally foraad dlnplad sheets and flat saptw sheets of tlvanluw foil. As da~ 
scribed in reference [3], this wultiwall concept laqMdas all three nodes of 
heat transfer— '—conduction, radiation and convection. The suparplastlcally 
fonad dlnplad sheets and tha long thin conduction path tend to alnlalsa heat 
conduction. Tha flat saptun sheets of tltaniua foil lapada radiation. The 
saall Individual voluaas created by tha dlaplad layers virtually allalnate air 
convection. Tha optlaal daslgn of such tharaal protactlon systaas requires 
effective techlquas In coupled tharaal and stress analyses. Finite alownt 
aethods offer tha greatest potential In aodallng such conpllcatad problaaa. 
However, tha resulting saal*-dlscrata aquations aay Invol^ asny thousand 
degrees of fraadoa. Since tha problem to be solved Is trsnsiant and ooo- 
llnaar, the selection of an appropriate tlaa Integration aathod la an essen- 
tial step In tha solution of such a coapllcatad problaa. Adalaan and Hafka 
[5] recently conducted a survey study on tha parforaanca of explicit and 
Implicit algorlthas for transient tharaal analysis of structures. Calcula- 
tions ware carried out using the SPAR finite <viaaant coaputar prograa [6] and 
a special purpose finite alaaant prograa Incorporating tha CBARB and (ZARIB 
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algorlchu. Bas«d upou th«ir studlMi thaj concluded thtt, gvmirallyt 
It algorithM mr« praf arable to explicit algorltlaw for "atiff problMM* 
though noir^convergence and/or «lde~bendlng of the resulting Mtrix equations 
■ay decrease the advantage of the Ispliclt ■etho’ls. 

These difficulties are aiailar to those found in fluid*~structure prob~ 
leas. Over the past few years, several reaedies have been proposed for these 
difficulties. Belytschfco and Mullen [7] have proposed an expllclt-iapliclt 
■ethod where the assh is partitioned into doaains by nodes and the partitions 
are sieultaneouly integrated by explicit and laplicit aethoda.' Hughes and 
Uu (8) have proposed an alternate iq>llcit*^xpllclt finite eleaent aethod 
where the aesh la partitioned Into d'.naina by eleaents and this eleaent parti*' 
tlon concept siiq>llfles the coi^>uter-iapleaentation and enhances its coapati- 
billty «flth the general purpose finite eleaent software. 

Although the laplicit*^xpllclt aetUod has been proven to be very auccesa- 
ful in soae fluid-structure Interaction probleas (see e.g., [8-10]), the site 
and coaplexlty of the prograa are increased because of the addition of the 
laplicit aethod. To overcoae these difficulties, Belytschko and Mullen [11] 
have proposed an partition, in which explicit tlae integration la used 
throughout. However, different tiae steps within different parts of the aesh 
can be eaployed slaultaneously. Partitioned and adaptive algorithas for ex- 
plicit tlaa integration have also been proposed by Belytschko [12]. 

Recently, Uu and Belytschko [13] put forward a general nixed tiae 
iaplicit-expllcit partition procedure within a linear context. It incorpo- 
rates the aentioned algorithas as special cases and is shown to have better 
stability properties than that in ^-S partition [11]. Siailar concepts esn 
also be used in transient conduction forced-convection analysis (see Uu and 
Un [14]). 
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In th« prcsnnt report, wo extend these lapllclt-explleit eoaeepte (nodes 
end elenente) to trenelent thenul enelysis of structures where different tine 
Integrstlon nethods with different tine steps esn be used In each elenent 
group* The sin of this approach la to achieve the attributes of the vsrlous 
tine integration nethods* 

For axanple. In transient structural analysis, explicit nethods require 
the size of the tine step to be proportional to the length of the shortest 
elenent; while In transient themal analysis, explicit nethods require the 
step size to be porportlonal to the the square of the length of the shortest 
elenent. So It Is wore advsiut«xcous to enploy this nixed tine Inpllclt-* 
explicit technique for transient themal analysis of structures since the 
partition proposed in [11,12] is often Inefficient for this kind of 
problen though It Is very efficient In structural analysis. 

In Appe?idix 1, we present the stability analysis of nixed tine 
integration schenes for transient themal analysis* This chapter is now 
published in Nunerical Heat Transfer Vol. 5, pp. 211-222, 1982. 

In Appendix 2, we present the co^iuter Inplenentation aspects of these 
nixed tine partition procedures. This chapter has been accepted for 
publication In International Journal of Nunerical Msthods In Engineering . 

In Appendix 3, we present the inprovenent of nixed tine Inpllclt-expllclt 
algorlthns for themal analysis of structures. In particular, three nunerical 
exanples are presented to evaluate the perfomance (l.e. accuracy and 
stability behavior, conputet storage and solution tine, etc.) of these nixed 
tim finite elenent algorlthns. This chapter has been accepted for 
publication in Conputer hethhods in Applied Mechanics and Engineering. 


In Appendix 4, we present a nethod for perfoming efficient and stable 
finite elenent calculations of heat conduction with quadrilaterals using one- 



point quadrature. Coaparlson with finite difference fonulaa haa shom that 
varloua values of the stabilised paraoeter« the S-polnt and 9**polnt aolecules 
can be obtained. It Is found that a combination of this one-point quadrature 
element and the mixed Itms Implicit-explicit methods may be an effective 
compromise. Ihls chapter has been submitted to International Journal of 
Numerical Methods In Engineering for possible publication. 

Interim report Is presented In chapter 5. 

In summary, the most Important mork In progress Is the development of 
nonlinear mixed time Integration methods for trenslent thermal analysis of 
stmctures suitable for Incorporation Into most finite element computer 
codes. Currently, continuous efforts are being made to Include three 
dimensional and nonlinear tbemal analysl.^. 

It should be emphaslxed that although the presently developed method has 
been applied to TPS (Thermal Protection Systems), the same technique could be 
applied to a vide variety of heat transfer fluid flow problems such as 
transient conduction forced-convection analysis, radiation and compressible 
fluid flov problems. 
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